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ABSTRACT 



AND Evan Scannapieco" 



During the early evolution of an AM CVn system, helium is accreted onto the surface of a white 
dwarf under conditions suitable for unstable thermonuclear ignition. The turbulent motions induced 
by the convective burning phase in the He envelope become strong enough to influence the propa- 
gation of burning fronts and may result in the onset of a detonation. Such an outcome would yield 
radioactive isotopes and a faint rapidly rising thermonuclear ".la" supernova. In this paper, we present 
hydrodynamic explosion models and observable outcomes of these He shell detonations for a range of 
initial core and envelope masses. The peak UVOIR bolometric luminosities range by a factor of 10 
(from 5 X 10^^ — 5 x 10^^ erg s~^), and the R-band peak varies from Mpt^peak = —15 to —18. The rise 
times in all bands are very rapid (< 10 d), but the decline rate is slower in the red than the blue due to 
a secondary near-IR brightening. The nucleosynthesis primarily yields heavy a-chain elements (^"^Ca 
through ^^Ni) and unburnt He. Thus, the spectra around peak light lack signs of intermediate mass 
elements and are dominated by Ca ii and Ti ii features, with the caveat that our radiative transfer 
code does not include the non-thermal effects necessary to produce He features. 

Subject headings: binaries: close — novae, cataclysmic variables — nuclear reactions, nucleosynthesis, 
abundances — supernovae: general — white dwarfs 



1. INTRODUCTION 

AM Canum Venaticorum systems (AM CVns) are 
He-transferring binaries with orbital periods of 5 — 60 
min consisting of a degenerate or semi-degenerate He 
don or and a C/0 or 0/Ne white dwarf (WD) accre- 
tor (jWarneii [19951 iNelemansI [20051) . There are ~ 20 
known sources, many discovered r ecent l y with the Sloan 
Digital Sky Survey (York ct al.' '2000'; Anderson et al.' 
[2005; Roelofs ct al. 2005, 2007, 2009; Ran ct al. 2010). 
There are sev eral proposed p r ogenitor scenarios for 
these systems (iPaczvhskil 119671: llben k, Tutukovl 119911: 
iPodsiadlowski et al.ll200l : however, recent work on the 
composition of the donors suggests that many, if not 
most, of the known AM CVns had degenerate He WD 
donors when mass transfer began (jRoelofs et al.l 120091 : 
[Nelemans et al.,,2009i) . 

In the He WD donor scenario. He mass transfer begins 
at rapid rates M = 1 — 3 x 10~^ Mq yr~^ and decreases 
with time as the binary evolve s under the influence 
of gr a vitational wave radiat i on (iTutukov fc YungelsonI 
T99a iNelemans et al.l [200lt iDelove fc BildstenI 120031 



Delove et al.l |2005|V During this initial phase of rapid 



accretion, the nuclear burning in the accreted He en- 
velope is thermally stable, and the binary may appear 
as a supersoft X-ray source with an orbital period of a 
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few minutes (Ivan Teeseling et"allll997l: IShen fc Bildst"CTl 
|2007[) . As M drops, the burning becomes thermally un- 
stable, yielding ~ 10 He flashes in the ^ 10^ yr before the 
amount of accreted mass required to achieve He ignition 
become s larger than the ava ilable fuel remaining in the 
donor ([Bildsten et al.ll2007D . Nearly aU observed AM 
CVns have evolved past this stage and are now slowly 
transferring He without any nuclear-powered phcnom- 
cna[3 

The largest ignition masses achieved during the He- 
flashing stage of AM CVn evolution in the He WD 
donor scenario are ^ 0.1 Mq, large enough to possi- 
bly yield hydrodynamical burning and He detonations, 
which would be potentially observable as faint ther- 
monuclear ".la" su pernovae (SNe .la; iBildsten et al.l 
120071: IShen fc Bildstc n 2009). In this paper, we fur- 
ther motivate our previous assumption of a He det- 
onation and calculate observables for these predicted 
events. In Sj2l we calculate the conditions at which 
the convective burning phase of the He flash becomes 
inefficient at transporting entropy throughout the en- 
velope. We find that the subsequent onset of hydro- 
dynamical burning occurs in the turbulence-dominated 
"d istributed" burning re g ime, which has been stiggeste d 
bv iNiemever fc WooslevI (|1997[ ). iKhokhlov et al.l |l997| ). 
and others to lead to a deflagration-to-detonation transi- 
tion in Type la supernovae. We thus proceed under the 
assumption that the He-burning yields a detonation. In 
121 we describe the hydrodynamic and nucleosynthetic 
evolution of these He detonations under the assumption 
of ID spherical symmetry. In §11 we use the Monte Carlo 

The short orbital period systems RX J0806, V407 Vul, and 
ES Cet may still be in this He-flashing stage, but their parameters 
are too uncertain to say som ething definit ive (Cropper ct al. 1993; 
IsraeLet al. "2002; War ner .STWoudg 12002: iStrohmaveril200a I2058l : 



Delove Taam 2006). 
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radiative transfer code SEDONA (jKasen et al.ll2006D to 
model the light curves and spectra from these events. We 
compare to known peculiar supernovae and conclude in 

m 

2. CONVECTIVE TO DYNAMIC PHASE 

During the accretion phase of the He flash cycle, radia- 
tive diffusion efficiently transports the entropy that leaks 
out of fluid elements within the accreted layer. However, 
when nuclear reactions at the base of the envelope be- 
come significant, radiative diffusion is unable to match 
the increase in entropy generation, and convection be- 
gins. The hydrostatic evoluti on of the growing convective 
zo ne ( Fuiimoto & Sugimotol 119821) . as recently outlined 
in iShen fc BildstenI (1200 9) , is long enough that spherical 
symmetry is a good assumption. 

The convective burning phase continues until convec- 
tion becomes inefficient at transporting the energy gen- 
erated by triple-a burning at the base of the envelope 
throughout the whole convective region. For an analo- 
gous analysis in Type la super no vae (SNe la), s e e, e.g 



Hillcbrandt k Nicn ieverl (|200 
Kuhlc n ct al. (200; 



Woosle 
and iPiro fc Chan 



et al. 



mm . 



This 

transition occurs when the local heating timescale at the 
base of the convective zone, theat = cpTt/eb, becomes 
shorter than the convective eddy turnover timescale, 
^eddy = H/vc, after which the evolution of the burning 
layer effectively decouples from the rest of the convective 
region. Here cp = SkTb/Snip is the specific heat at con- 
stant pressure in electron degenerate He, 6-subscripts re- 
fer to the envelope base, Cb is the energy generation rate, 
H = Pb/pbg is the pressure scale height, and g is the local 
gravitational acceleration. In mixing-length theory, the 
convective velocity, Vc, can be expressed via the convec- 
tive flux as Fc £bPbH ^ PbVc^ so that Vc ~ {ebHy^^. 
When evaluated at ihcat = ^cddy, the convective ve- 
locity in electron degenerate He is v'^ ^ {cpTbY^"^ '-^ 
2 X 10** cm s Tg , where Tg is the base temperature 
in units of 10^ K at the onset of dynamic burning. This 
velocity is roughly equivalent to the thermal speed of the 
ions, but both are still subsonic because the sound speed 
is determined by the degenerate electrons and is a factor 
of (Ep /kXy^^ faster, where Ep is the Fermi energy; see 
column 4 of Table [5] for the initial degree of degeneracy 
for our models. 

This large scale convective motion yields a cascade of 
turbulent velocities to shorter lengthscales, I; if we as- 
sume Kolmogorov scaling, these velocities are given by 
fturb(0 ~ Vc{l/ HY^^ ■ After the transition to inefficient 
convection, the bulk convective motion and its associated 
turbulent cascade is effectively frozen out with respect to 
the evolution of the burning layers. Thus, there is a spec- 
trum of turbulent velocities present at the initiation of 
any burning wave in the convective shell, independent 
of the buoyancy-driven turbulence associated with rising 
flames 

* That said, the convcctively-driven turbulent cascade does not 
differ strongly from the cascade driven by a buoyant flame. The 
buoyancy speed of a laminar He flame is shown in Figure 2 of 
ITimmes fc Niemevej JJOOO) to be 3 X lO'' - 10* cm s"\ which is 
roughly equal to = 8 X 10'' cm s^ ^ , calculated at T;, = 2 X 10® K, 
and since the lengthscale for both processes is ~ the turbulent 
cascades are nearly the same. 



The presence of these turbulent motions will dramat- 
ically modify the propagation of a laminar deflagration 
wave if the turbulent velocity at the lengthscale of the 
laminar flame width exceeds the burning speed. Presum- 
ing just the turbulence generated by buoyancy (but see 
footnote 1), studies of laminar li e deflagrations ( Timniei 
120001: iTimmes fc NiemeveHl2000t) have found that turbu- 
lence indeed dominates the flame properties for initial 
P ^ lO'' g cm^^. This turbulence-dominated condition 
is often referred to as the "d istribut ed" burning regime 
fepc 1987- Niemever fc Wo oslev 1997t iPetersI I2OOOI : 
lAspden et all 120081 ). AM CVn evolution yields convec- 
tive burning envelopes with densities o f < 3 x 10^ g cm~'^ 
(|Bildsten et al.l2007t[Shen fc Bildstel3l2009| ) and the pro- 
longed convective phase ensures a turbulent medium, so 
that deflagrations in these He shells will be in the dis- 
tributed burning regime. 

In the SN la literature, it has been suggested 
that a deflagration in the distributed regime can 
subsequently transition to a detonation via the 
Zel'dovich gradient mechanism ([Z el'dovich et al.l 11970 ; 



Blinnikov fc Khokhlov 1986, 1987; Khokhlov ct al. 1997 



Niemever fc Kcrstein 1997; Nie mever fc W ooslev ,1997 



Seitenzahl et al.ll2009f ). One criterion for the transition 



to a detonation is that the large scale t urbulent veloci- 
ties should not be to o subsonic (> 0.2cc: IWooslevI 120071 : 
iWooslev et al.l [2009t ) . For a convectively- driven turbu- 
lent spectrum fixed approximately at the point that con- 
vection becomes inefficient, the convective velocity is 
v'^ ~ 0.3cs (a factor of yJkT/Ep slower than the sound 
speed as discussed previously), which meets this crite- 
rion. The exact outcome of distributed burning in SNe la 
is still unclear (Niemeycr 1999), and the conclusions are 
even muddier in the case of less well-studied He defiagra- 
tions, but it is plausible that a detonation will result from 
convectively-inefficient burning in AM CVn systems; for 
the rest of this study, we will examine the outcome of 
such shock-driven burning, with the caveat that a He 
deflagration may yield a qualitatively different outcome. 

We will further assume a spherically symmetric ge- 
ometry for the system even though multi-dimensional 
effects will be important for the explosion. In partic- 
ular, the detonation will likely be initiated at a single 
point and spread tangentially around t he WD with quasi- 
cylin drical symmetry, as explored by iFink et aTl (|2007l 
|2010| ). We leave this complication for future work. 

3. HYDRODYNAMIC EVOLUTION AND 
NUCLEOSYNTHESIS 



In iBildsten etall (|2007l ) and iShen fc BildstenI (|2009l) . 
we performed preliminary calculations for the initial con- 
ditions and general appearance of these He shell detona- 
tions. We now proceed to probe the range of observable 
SN .la outcomes. We evolve six combinations of core 
and envelope masses (core -I- envelopes masses: 0.6-1-0.2, 
0.6-h0.3, 1.0+0.05, 1.0+0.1, 1.2+0.02, and 1.2+0.05 Mg) 
from the onset of hydrodynamical burning. Note that en- 
velope masses > 0.1 Mq are not expected to be achieved 
during the evolution of AM CVn systems in the He WD 
donor scenario; however, we include them in our mod- 
els to allow for their existence in the He-burning star 
donor scenario (Iben fc Tutukwl[T99lt see, e.g.. Fig. 1 
of ,Shen fc Bildsten .20091) . The lower of the two enve- 
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Fig. 1. — Snapshots of the temperature profile (top panel) and 
pressure profile {bottom panel) of the 0.6 + 0.2 Mq model during 
its hydrodynamic evolution. Solid lines show the temperature and 
pressure vs. the enclosed mass at 0, 0.08, 0.2, 0.3, 0.4, and 0.5 
s after the beginning of the hydrodynamic simulation when the 
temperature at the base of the envelope is perturbed upward by 
~ 10%. 

lope masses for each core mass is chosen to be near 
the minim um enve lope mas s for hyd rodynamical burn- 
ing (Bil dsten et al.|[2007: Shen fc Bilcisten.200 9'). and the 
higher envelope's mass is chosen to be roughly twice as 
large. 

The isothermal T — 10^ K cores are composed of 50% 
^^C and 50% ^^O by mass, and the initially hydrostatic 
and isentropic envelopes are composed of pure He. The 
initial thermodynamic conditions at the envelope base 
are set near the point when burning becomes hydro- 
dynamic, as described in [J2] This initial configuration 
is evolved with a ID spherically symmetric explicit La- 
grangian hy drodynamic code (see, e.g. , [Bcuz 1991]), as 
described in lWeinberg fc BildstenI (|2007f ) . The 13 isotope 
a-chain nuclear networlfl includes a-capture, heavy-ion, 
and (a,p)(p,a) reactio ns. Different initial compositions 
(jShen fc BildstenI |2009() and neglected nuclear reactions 
could change the nucleosynthetic outcome of these mod- 
els; we leave an exploration of these effects for future 
work. 

At the beginning of the hydrodynamic simulation, the 
temperature at the base of the envelope is perturbed up- 
wards by ^ 10%, which results in a thermonuclear run- 
away and a pressure perturbation that steepens into a 
shock and an outwardly propagating detonation. Figure 
[T] shows the evolution of the temperature and pressure 
profiles for the 0.6 -I- 0.2 M© model during this phase. 
The solid lines mark the temperature and pressure vs. the 
enclosed mass (i.e., the surface is to the right) at 0, 0.08, 
0.2, 0.3, 0.4, and 0.5 s after the beginning of the hydro- 
dynamic simulation. 

After the shock reaches the surface of the envelope, 
nuclear reactions are turned off for computational sim- 
plicity (as the post-shock breakout nucleosynthesis was 

^ http: / / www .cococubed.com/code_pages / net_aproxl3.shtml 
and ri'immea ||T999| ') 
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Fig. 2. — Free-streaming ejecta velocity vs. mass above for all 
models: 0.6 + 0.2, 0.6-1-0.3, 1.0-1-0.05, 1.0-1-0.1, 1.2-1-0.02, and 
1.2 + 0.05 Mq. Core masses are as labeled; black, red, and blue 
lines are for core masses of 0.6, 1.0, and 1.2 Mq, respectively. The 
solid lines are the smaller envelope masses for each model, and the 
dashed lines are the larger envelope masses. The velocities shown 
are calculated by taking the velocities 100 s after shock breakout 
and subtracting the effect of leaving the gravitational potential 
wells. 



TABLE 1 
Ejecta parameters 



MwD 

[Mq] 


Monv 

[Mq] 


[Mq] 


^kin,oo 

[10^" erg] 


«>^^^ 
[lO"' km s-i] 


^dur 

[d] 


0.6 


0.2 


0.20 


2.7 


1.2 


6.6 


0.6 


0.3 


0.30 


5.3 


1.3 


7.6 


1.0 


0.05 


0.048 


0.53 


1.1 


3.4 


1.0 


0.1 


0.090 


1.2 


1.2 


4.4 


1.2 


0.02 


0.017 


0.13 


0.90 


2.2 


1.2 


0.05 


0.036 


0.43 


1.1 


2.9 



found to be insignif icant in a test case run with FLASlf^: 
iFrvxell et all [2000l ). The core and envelope are fur- 
ther evolved until the ejected material reaches a free- 
streaming, homologously expanding state. Figure [5] 
shows the free-streaming ejecta velocity vs. mass above 
(i.e., Mabovc = Mtotai " Mc„c) for all models: 0.6 + 0.2, 
0.6 + 0.3, 1.0 -I- 0.05, 1.0 -I- 0.1, 1.2 + 0.02, and 1.2 -|- 
0.05 Mq. The velocities shown are the velocities 100 
s after shock breakout with a small correction to ac- 
count for leaving the gravitational potential well: i.e., 
v^- = — 2GMcnc/r, where the velocity, v, and radius, 
r, of the fluid elements are evaluated at 100 s after explo- 
sion. We neglect any impact of the binary motion. For 
larger core masses and thus larger initial gravitational 
potential wells, less of the envelope is able to escape, so 
that nearly 100% of the 0.2 M© envelope on the 0.6 M© 
core is ejected, whereas only ~ 75% of the 0.05 M© en- 
velope on the 1.2 M© core is ejected. 

The ejecta parameters for each model are summa- 
rized in Table [U which shows the ejecta mass, M^j (col- 
umn 3), total integrated kinetic energy in the ejecta, 

http;//flash. uchicago.edu 
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kin,! 



2 \l/2 = 



^ {v1-j2)dm (column 4), average ejecta veloc- 
■^2£'kin,oo A'^cj (column 5), and light curve 



ity, («cj) 

duration, idur = \/ KMcj/7c(uoj), with assumed opacity 
K = 0.2 cm'^ g""'^ for a He-rich medium (column 6; see 
g]for more details). The average ejecta velocities are all 
very similar, between 0.9 — 1.3 x 10* km s~^. This re- 
sults from the similar mass fractions of burned material 
in each model as well as the fact that burning from *He 
to '^^S releases only 10% less energy per mass as burning 



all the way to ^^^Ni (1.4 x 10^** vs. 1.5 x 10^** erg g-^); 
see !j3]and Table [2] for more details. 

As evident in Figure [U the strong pressure in- 
crease due to the hydrodynamic burning in the en- 
velope sends a shock wave into the WD core. A 
sequence of events like this has been discussed ex- 
tensively in the context of double degenerate Type 
la supernova progenitors, either as an edge-lit deto- 
nation in which the s hock ignites C-burning at the 
core-envelope interface (iNomotol 11982 1: iLivn e fc Glasner 
1990, 1991; Livne & A rnettl ll995n Wiggins et all I199S ' 
Garcia- Senz et al.i il999). or as a shock- focused detona- 
tion in which the geometrically-focused shock yi elds a 
carbon detonation near the center of the WD (iLivnd 
19901: iLi^ fc GlasneH [T99lt IWooslev fc WeavedTT99r 



Livne fc Arnettl 1199,1 iFink et all 120071 120101 : ISim eTaLl 
20101 ). In our models, where the temperature is only 
perturbed at the exact interface between the envelope 
and core, the inwardly propagating shock is not strong 
enough to immediately ignite the ^^C at the core's 
edge. For detonations initiated well above the base of 
the envelope, the inwardly propagating shock strength- 
ens as it approaches the envel ope-core interface and 
may ignite the edge of the core ([Livne fc Glasnedll990l: 
^arcia-Senz et al. 1999); we do not consider this possibil- 
ity. Given the uncertainties in these detonation models, 
and the expectation that some fraction of the AM CVn 
population contains 0/Ne WD accretors, which would be 
much more difficult to detonate via these mechanisms, 
we proceed under the assumption that the core is not 
detonated, with the caveat that this outcome remains a 
possibility. 

Table [2] shows the ejecta mass (column 5) for each 
model, denoted by its core and envelope masses (columns 
1 & 2), initial base density (column 3), and initial degree 
of degeneracy (column 4), along with the mass fractions 
of ^He, 36 Ar, "OCa, '''*Ti, "SCr, ^Spe, and in the 

ejecta (columns 6-12). Column 13 shows peq, the density 
where the He and Ni mass fractions are equal, except for 
the 1.0-1-0.05 Mq model for which the two mass fractions 
were not equal anywhere in the envelope. For the larger 
envelope masses for each core mass, the majority of burn- 
ing products in the He shell are in the form of ^^Ni, with 
decreasing contributions from ^^Fe, ''^Cr, and ^''Ti. For 
the smaller envelope masses, the contribution of isotopes 
below ^^Ni are enhanced as compared to the larger enve- 
lope masses. In particular, the 1.0+0.05 Mq model, with 
a relatively low initial base density ~ 6 x 10^ g cm"^, 
only produced a small amount of ^^Ni. Its nucleosynthe- 
sis is dominated by the production of '^''Ca. 

Profiles of mass fraction vs. ejecta mass coordinate 
are shown in Figure [3] for a representative model (0.6 + 
0.2 Mq). Solid lines are for the isotopes as labeled, and 




Fig. 3. — Mass fractions vs. enclosed mass at the time of shock 
breakout for the 0.6-1-0.2 Mq model. Solid lines are for the isotopes 
as labeled, and dashed lines are for the remaining isotopes, which 
are insignificant in the ejecta. 



dashed lines are for the remaining isotopes, which are 
insignificant in the ejecta. Near the base of the initial 
He shell, where the fluid elements are confined at high 
temperatures and densities for the longest time, the nu- 
cleosynthesis is dominated by Fe-group elements at the 
end of the a-chain. Closer to the surface, where there 
is little time to burn the fuel prior to radial expansion, 
most of the He remains unburnt. For the densities and 
temperature reached in these detonations, the slowest 
step of He-burning is the first: triple-a burning. If there 
is time to form a ^^C nucleus, it will capture additional 
He nuclei relatively rapidly and proceed all the way to 
Fe-group elements near the end of the a-chain. Thus, for 
all the models, very few intermediate mass elements are 
produced; the cumulative mass fraction of isotopes be- 
tween ^°Ne and ^^S for each model is < 1%. iPerets et"aLl 
(|2009[ )'s one-zone analysis produces a qualitatively sim- 
ilar nucleosynthetic yield for the relevant temperatures 
and their most He-rich composition. This has strong 
implications for the spectral signatures of these models, 
which will be discussed in !i4.2l 



4. RADIATIVE TRANSFER AND OBSERVABLES 

After reaching the free-streaming phase, the hydrody- 
namic output described in ij3]is process ed by SEDO NA, a 
Monte Carlo radiative transfer code fiKasen et al.l l2006). 
The necessary modifications to follow the radioactive de- 
cays of ""^^Cr and ^^F e have been included as outlined in 
iBildsten et all pOOTi ). 

We can gain a qualitative estimate of the light curve 
properties by employing the simplifying assumptions de- 
rived in lArnetti (l982) and Pinto fc Eastma n (2000a): 
namely, that the duration of the light curve, tdur, has 
a scale set primarily by the time at which the radiative 
diffusion time equals the age of the explosion (column 
6 of Table [1]) , and that the peak bolometric luminosity, 
ipeak, is equal to the instantaneous radioactive decay at 
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TABLE 2 

EjECTA COMPOSITION 



MwD 

[Mq] 


[Mq] 


Pb 




[Mgj] 


Xi 


^36 




X44 


^48 


X52 




Poq 

[10^ g cm"-'] 


0.6 


0.2 


6.8 


0.28 


0.20 


0.45 


0.0067 


0.036 


0.12 


0.13 


0.12 


0.13 


4.1 


0.6 


0.3 


24 


0.083 


0.30 


0.30 


1.2 X 10"" 


0.0020 


0.0072 


0.0097 


0.022 


0.65 


3.6 


1.0 


0.05 


5.8 


0.40 


0.048 


0.44 


0.034 


0.30 


0.11 


0.032 


0.033 


0.015 




1.0 


0.1 


21 


0.10 


0.090 


0.39 


5.8 X 10~5 


0.0014 


0.0066 


0.011 


0.032 


0.56 


5.1 


1.2 


0.02 


9.3 


0.28 


0.017 


0.51 


0.0037 


0.025 


0.063 


0.081 


0.17 


0.13 


5.8 


1.2 


0.05 


34 


0.068 


0.036 


0.38 


3.6 X 10"'' 


8.0 X 10"" 


0.0033 


0.0052 


0.015 


0.59 


6.6 



t = tdurEB These basic assumptions yield 



and 



0.2 cm^ g- 



0.1 Mq 



10 km s 



1/2 
(1) 



-'peak 



■ exp 



n 



(2) 



where the i-subscript refers to the dominant radioactive 
product in the ejecta of which there are Ni nuclei, Qi is 
the energy release per radioactive decay, and is the e- 
folding lifetime for the decay. This approximation of the 
peak luminosity yields Lpcak ~ 2 x 10"^^ erg s"^ or a peak 
bolometric magnitude Mboi ^ -17 for 0.05 M© of ^'^Ni 
and tdm = 5 d. Note that this estimate assumes that the 
ejecta is optically thick to all radiation, including 7-rays, 
at peak. This estimate is further complicated by the 
inclusion of multiple radioactive burning products, each 
with more than one step in their decays to their stable 
end states. Additionally, these radioactive isotopes have 
different diffusion times due to their vertical stratification 
in the envelope, as shown in Figure[3] However, equation 
([2]) is useful in describing general trends in Lpcak- In 
particular, it demonstrates two offsetting features of SNe 
.la: while the low radioactive mass (i.e., small Ni) lowers 
-^pcak as compared to SNe la, the low ejecta mass implies 
a short tdur, thus increasing Lpcak- The net effect is to 
decrease Lpoak, but not as much as one might naively 
assume given the amount of radioactive material. 

4.1. Light curves 

Figure |4] shows UVOIR bolometric light curves for the 
six models as labeled. Solid lines are for the lower en- 
velope mass for each core, and dashed lines are for the 
higher envelope mass. Black, red, and blue lines are for 
core masses of 0.6, 1.0, and 1.2 Mq, respectively. Squares 
denote the time after which over half of the gamma-rays 
escape the ejecta without bei ng absorbed. As exp ected 
from the above discussion and lBildsten et al.l (|2007f ). the 
light curves reach peak luminosity fairly rapidly, within 
2 — 10 d after the detonation, with peak luminosities in 
accord with the energy input from instantaneous radioac- 
tive decays. The peak luminosities range over an order of 
magnitude from 5 x lO"'^ — 5 x 10**^ erg s~^, corresponding 
to peak bolometric magnitudes of —15.5 to —18. 

Radioactive energy deposited in mass shells closer to 
the surface have shorter diffusion timescales. Thus, the 

The low mass envelopes remain radiation dominated through- 
out their evolution. 



,0.6 + 0.3 




5 10 15 

Time since explosion [d] 

Fig. 4. — UVOIR bolometric luminosity vs. time since explosion 
for all models as labeled. Solid lines are for the lower envelope 
mass for each core, and dashed lines are for the higher envelope 
mass. Black, red, and blue lines are for core masses of 0.6, 1.0, and 
1.2 Mq, respectively. The multiple peaks are due to the stratifica- 
tion of radioactive isotopes in the ejecta. Squares denote the time 
after which over half of the gamma-rays escape the ejecta without 
being absorbed. 



stratification of radioactive isotopes in the ejecta, as 
demonstrated in Figure [3] for the 0.6 -|- 0.2 model, yields 
multiple peaks in the bolometric light curves of models 
with envelope masses > 0.1 Mq for which the light curves 
are broad enough to distinguish the multiple peaks. In 
these models, the first peak or shoulder is due to the 
rapid decay of ''^Cr and ^^Fe near the surface, followed 
by a broader peak due to the decay of these isotopes and 
^^Ni deeper within the ejecta. 

The Vega magnitudes in the UBVRI filters vs. time 
since explosion for the 0.6 -I- 0.2 Mq model are shown 
as solid lines as labeled in Figure (5] the magnitude of 
the bolometric light curves is shown as a dashed line. 
The early peak in the bolometric light curve due to the 
stratification of isotopes ~ 2 d after the explosion yields 
an early peak in the U filter, since the luminosity, radius, 
and the Wien displacement law give a maximum in the 
spectral energy at ~ 2000 A at this time. The second 
bolometric peak, primarily powered by the decay of ^^Ni, 
yields maxima in the other filters at roughly the same 
time as each other, ~ 8 d after explosion. 

Figure [6] shows the R-band light curves {top panel) and 
V — R colors (bottom panel) for all models. Bullets mark 
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Fig. 5. — Multiband light curves for the 0.6 + 0.2 M© model 
in Vega magnitudes. The dashed line shows the bolometric light 
curve, and the solid lines show magnitudes for various filters as 
labeled. The secondary peak in redder wavebands is apparent be- 
ginning ~ 15 d after explosion, which is ~ 7 d after the time of 
B-band maximum. 
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Fig. 6. — R-band light curves {top panel) and V—R colors {bottom 
panel) for all models as labeled. Solid lines show light curves and 
colors for the smaller envelope for each core mass, and dashed lines 
represent the larger envelopes. Black, red, and blue lines are for 
core masses of 0.6, 1.0, and 1.2 Mq, respectively. Bullets denote 
the time of the B-band maximum for each model. 



the time of the B-band maximum for each model. All 
of the models have V — R = —0.2 to at the time of 
B-band maximum, and all then become much redder at 
later times, with most possessing a secondary peak in 
the R, I, and redder bands. For the 0.6 -I- 0.2 Mq model 
shown in Figure [51 for example, the B-band light curve 
takes 6 d to drop 1 mag from B-band peak, while it takes 
12 d for the R-band light curve to drop 1 mag from its 
peak. This strong reddening is due to the same physics 
that produces secondary near-IR peaks in typical SNe la: 
namely, the emissivity of Fe-group isotopes peaks sharply 
in near-IR filters when they transition from the doubly- 
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Fig. 7.— Optical spectra for the 0.6-1-0.2 Mq model at -5, +0, 
-1-5, and -|-10 d with respect to B-band maximum, which occurs 8 d 
after explosion. The spectra are dominated by Ti ll features, along 
with the Ca ll H&K and IR lines and a short-lived Cr II feature 
near 4600 A. The lack of Si II features is also a notable signature. 
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4.2. Spectra 

The optical spectra of the 0.6 -I- 0.2 Mq model are 
shown in Figure [7] at —5, -1-0, -1-5, and -1-10 d with re- 
spect to the B-band maximum, which occurs 8 d after the 
explosion. The spectra are offset by arbitrary constants. 
Line and feature identifications in the spectra are made 
by comparing the complete spectra to radiative trans- 
fer simulations in which the effect of a single isotope is 
removed. 

The pre-maximum spectrum is a fairly featureless 
blue continuum. By the time of B-band maximum, Ti 
II features have become prominent. The Ca ii H&K 
lines and IR triplet are also present with blueshifts of 
1.4 X 10^ cm s~^, as is a Cr ii feature near 4700 A that is 
likely unique to supernovae powered by He-burning prod- 
ucts. At 10 d after the time of B-band maximum (18 d 
after explosion) , the continuum has become much redder 
due to the reprocessing of UV and blue light by Fe-group 
elements, as previously discussed. The Ca ii lines are still 
visible, but the spectrum is determined primarily by Ti 
II features. The Cr feature at 4700 A has disappeared by 
this point due to ^^Cr's short 22 hr half-life. 

The spectra are also notable in what they do not show. 
First, they lack lines from intermediate mass elements 
such as Si, S, and Mg. This is not surprising given the 
negligible amount of these isotopes produced in these 
models. Second, the spectra show no signs of He, which 
may seem counterintuitive given the large amount of un- 
burnt He in the ejecta. However, the relevant optical 
He lines have lower levels with high excitation energies 
(> 10 eV) and will not be highly populated by thermal 
processes given the moderate ejecta temperatures ~ 10* 
K. As the present calculations assume local thermody- 
namic equilibrium (LTE) to calculate the level popula- 
tions, none of the optical He lines has significant opac- 
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Fig. 8. — Optical spectra at B-band maximum for all models as 
labeled, with arbitrary offsets. Black, red, and blue lines are for 
0.6, 1.0, and 1.2 Mq cores, respectively. Solid lines are for the 
lower envelope mass for each model, and dashed lines are for the 
higher envelope mass. Times shown are the ages since explosion. 

ity. However, the gamma-rays from radioactive decay 
will Compton scatter and produce fast electrons that 
may non-thermally excite these levels and thereby in- 
crease the line opacity. Inclusion of this non-LTE effect 
is necessary to explain the He lines observed in Type lb 
supernovae (jLucy 1991). 

The strengths of the He lines in supernovae are there- 
fore not only sensitive to the He mass, but also to the 
amount of radioactive material present and on how this 
material is mixed into the He layer. In the SN .la case, 
radioactive material is produced adjacent to the He and 
the densities are relatively low, so gamma-rays do effec- 
tively diffuse into the He layer. It therefore seems likely 
that non-thermal excitation will produce He lines in some 
models, but we need quantitative non-LTE calculations 
in the future to confirm this. 

Figure |8] shows optical spectra for all the models as 
labeled at the time of their respective B-band maxima. 
Black, red, and blue lines are for 0.6, 1.0, and 1.2 AIq 
cores, respectively; solid lines are for the lower envelope 
mass for each model, and dashed lines are for the higher 
envelope mass. The times shown are the ages since ex- 
plosion for each model at their B-band maxima. Most 
models show Ca ii and Ti ii absorption lines with vary- 
ing strengths depending on the nucleosynthetic output 
and the photospheric temperature, which is itself a func- 
tion of the luminosity and radius at B-band maximum. 
E.g., the lower envelope mass models for each core mass 
(0.6 + 0.2, 1.0 -I- 0.05, and 1.2 -I- 0.02 Mq) are the most 
"^^Ca- and ■*'*Ti-rich, and two of these three models reflect 
this; however, the B-band maximum for the 1.0-1-0.05 Mq 
model occurs so early that the photosphere is hot enough 
to produce a relatively featureless blue continuum. 

5. CONCLUSIONS 

Previous wor k by iBildsten et al.l ()2007D and 
iShen fc BildstenI ()2009D hypothesized that a new 
class of SNe .la powered by He shell detonations could 



be produced by AM CVn systems, and in this work, 
we have calculated their observable properties in depth 
for the first time. In fj2l we examined the outcome of 
turbulent combustion with guidance from work done on 
Type la supernovae and found that the assumption of a 
He detonation is well-motivated by the initial conditions. 
The hydrodynamic evolution of these detonations was 
calculated in ID and led to the ejection of most of the 
envelopes with velocities of ^ 10'' km s~^ (Figs. [l]and 
[2] and Table [T]). The nucleosynthetic output of each 
envelope was dominated by unburnt He and a-chain 
elements from ^"Ca to ^^Ni (Fig. Hand Tabled]). We 
employed a Monte Carlo radiative transfer code to 
produce light curves (Figs. [H [5l and [6]), which are 
notable for their fast rises of 2 — 10 d, their multiple 
peaks due to the stratification of radioactive isotopes, 
and their strong late-time reddening due to Fe-group 
ionization effects. The spectra are dominated by Ca ii 
and Ti ii features (Figs. El and [5]). 

Recently, several peculiar Type I supernovae have been 
observed that do not fit cleanly into existing supernova 
categories. SN 2002bj was discovered in archiva l Lick 
Obser vatory Supernova Search data by Poznansk i et al.l 
(j2010( ). Its most striking feature is its extremely rapid 
rise and fall, with a > 3.5 magnitude R-band rise in 7 
d to its peak at ~ —18.5 in all bands, and a subsequent 
decline of 3 magnitudes in 15 d. A spectrum taken 7 
d after discovery shows signatures of He and possibly 
even V, which could be the decay product of ^^Cr. How- 
ever, while the general speed of the light curve and the 
evidence of He-burning are in line with the SN .la pre- 
dictions in this paper, SN 2002bj also shows evidence of 
intermediate mass elements and does not have Ti fea- 
tures or strong late-time reddening. 

SN 2008ha was a faint (peak My = -14.2) and fast 
rising and declinin g supernova (|Folev et al.l 120091 IMfll : 
iValenti et al.l 120091 ) that at first glance seemed a possi- 
ble SN .la. However, its slow ejecta velocity of ^ 2000 
km and the nucleosynthetic signatures of C-burning 
point to the failed deflagration of a C/0 WD instead of 
a He shell detonation. 

SN 2005E was a faint (peak Mb = -14.8) Ca-rich Type 
lb supern ova with a rapid r ise and a fall of 1 mag per 
10 — 15 d (jPerets et al.ll2009l) . Many of its characteristics 
match our general results, including a lack of interme- 
diate mass elements, although its bolometric light curve 
may rise too slowly (Fig. The nucleosynthesis does 
not match exactly, as SN 2005E shows evidence for a 
large amount of C a as well as some C and O, but given 
the nucleosynthetic outcome of our weakest explosion 
model (1.0-)-0.05 Mq), it is possible that a still weaker ex- 
plosion would yield an ejecta composition matching that 
of SN 2005E. As previously mentioned, the inclusion of a 
larger nuclear network and multidimensional effects will 
likely play a role in changing the nucleosynthetic yield; 
these refinements await future work. 

Figure [9] compares our model light curves to those of 
several observed SNe. UVOIR bolometric light curves 
are plotted as in FigurelU but vs. the time since bolomet- 
ric peak instead of the time since explosion. Thick lines 
show .la models as throughout this paper: thick solid 
lines show the less massive envelopes for each core mass, 
thick dashed lines show the more massive envelopes, 
and black, red, and blue lines denote 0.6, 1.0, and 



8 



SHEN, KASEN, WEINBERG, BILDSTEN, & SCANNAPIECO 



o 



o 











y^fiA (la) 


v02bj 




y' 






/ 99b>- 


Xx ■ 














r 














/ -m 


^ ■ 

1,0 •! d"""*..^ 






1-2 V 0.05 




■' I'' 1 ' W 




■ I 


' \' ' \ \ 

> !■ ( 1 \ 

1 li i-^'Tn8h^~---^ 

' M 1 

< 1 • / ' 1 

. i. 1 • . i/. 1 . J 1 . . . . 


\l2 + 0.02 ^'-l 



-10 —5 5 10 

Time since bolometric peak [d] 

Fig. 9. — UVOIR bolometric luminosity vs. time since bolometric 
peak. Our .la models are shown as thick lines as labeled; thick solid 
lines show the less massive envelopes for each core mass, thick 
dashed lines show the more massive envelopes, and black, red, and 
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are quasi-bolometric light curves for the three peculiar 
SNe just disc ussed as well as S N 2005cf, a well-studied 
bright SN la (IWaiig et all 120091) . and SN 1999bY, a sub- 
luminous SN la (jGarnavich et aD 120041 : iPhillips et "aD 
[2OOI). The light curves are constructed from integrated 
and blackbody fits to broadband photometry; in the case 
of SN 2005E, which was limited to BVRI observations, 
the U-band and near-IR contributions are assumed to 
scale as they do in SN 2005cf (S. Valenti, private com- 
munication; IPastorello et~all 12007) Fl 

It is clear that the predicted SN .la rise times of 2 — 10 d 
are much faster than any other rise times, aside from the 
extremely fast SN 2002bj, but the peak bolometric lumi- 
nosities cover a large range from those of faint peculiar 
SNe to sub-luminous SNe la, and the declines in lumi- 
nosity are not strongly atypical of observed SNe. Thus, 
the most robust and unique signatures of a SN .la will be 
a very fast rise to peak luminosity and the appearance 
of He-burning products such as Ca 11 and Ti 11 in the 
spectra. 

We thank Ryan Foley, Hagai Perets, Anthony Piro, 
Dovi Poznanski, and Dean Townsley for helpful discus- 
sions, and Stefano Valenti for providing the bolometric 
light curve of SN 2005E. This work was supported by the 
National Science Foundation under grants PHY 05-51164 
and AST 07-07633. 
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